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Abstract 

We have looked into an experiment that has been termed the "canonical example" of jamming: 
granular material, clogging the outlet of a container as it is discharged by gravity. We present 
quantitative data of such an experiment. The experimental control parameter is the ratio between 
the radius of the orifice and the radius of the beads. As this parameter is increased, the jamming 
probability decreases. However, in the range of parameters explored, no evidence of criticality 
— in the sense of a jamming probability that becomes infinitely small for a finite radius — has been 
found. We draw instead a comparison with a simple model that captures the main features of the 
phenomenon, namely, percolation in one dimension. The model gives indeed a phase transition, 
albeit a special one. 

PACS numbers: 45.70.Ht, 45.70.Mg 
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The suggestive idea that jamming is at the origin of a new state of matter, amenable to be 



described in a thermodynamic formalism, has been recently proposed Our experiment 
is akin to the familiar experience of shaking a saltcellar; salt grains plug the outlet of the 
container due to the formation of arches. Shaking or tapping is needed to break up the 
arch and restore the flow after each jam. A similar situation is found in many industrial 
applications, such as dosage hoppers, and in various transport phenomena. Ad hoc laws 
have been put forward based on experience and observation 2], but the physics of arching 
is still poorly understood. A new vision of jamming has recently been put forward, relying 
on the concept of fragile matter, which encompasses not only granular materials but foams, 
unstructured glasses and other systems The common feature they share is that they 
support certain stresses, called compatible, developed as a response to the external forces, 
but are unstable against incompatible stresses. 

We have devised a simple arrangement to conduct a quantitative study of the jamming 
that halts a granular flow. Our assembly consists of a scaled cylindrical silo with a circular 
opening at the base. The bin is filled with monodisperse glass spheres (radius standard 
deviation is about 1%). Beads of different sizes were used in separate trials, but always with 
a diameter larger than 1 mm in order to reduce the relative importance of disturbances - 
such as humidity or electrostatic interactions - against gravity. The diameter of the bin is 
larger than 30 bead diameters; in this way, the finite size of the system can be neglected, as 
reported in the literature J] and checked specifically for our set-up j^. The silo is always 
filled with the same procedure (a distributed filling 6], pouring the grains rapidly in the 
silo) in order to obtain a repeatable compaction of the material (a volume fraction of about 
0.59) from run to run. Besides, the pressure at the base of a silo is known to be almost 

Q 

constant as long as the height of the material above it exceeds a certain level [7(. Indeed, 
the results of our experiment are repeatable -irrespective of the amount of material in the 
silo- provided that the height reached by the beads is larger than twice the diameter of the 
silo. We refill the silo from time to time so that this condition is met. We have used a 
transparent glass silo in order to check that a mass flow develops inside it Q| ; a funnel was 
therefore not formed during the experiments. 

When the experiment begins, beads pour freely from the outlet due to gravity. Soon 
thereafter, however, if the diameter of the orifice does not exceed a few bead diameters, 
an arch is formed arresting the flow. An electronic balance beneath the silo measures the 
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number of grains s which fall between two successive jams. The size of the avalanche is 
stored in a computer and then the arch at the outlet is destroyed so as to trigger another 
avalanche. This is accomplished by means of a jet of compressed air aimed at the orifice from 
beneath the silo. In this way, the compaction of the granular material remains approximately 
constant during an experimental run - shaking would instead pack the grains together more 
densely |8|. Further information on the set-up will be provided elsewhere. We made sure 
that variables such as the pressure of compressed air or the material from which the bin is 
made have a negligible effect in the data that we are presenting [9|. The parameter that 
we control is the ratio between the diameter of the outlet orifice and the diameter of the 
beads, called R in the following. Both the size of the orifice and the size of the beads may 
be changed; only the ratio R is relevant. 

For a fixed value of R, the data collected after a large number of avalanches - several 
thousands typically - are best viewed in the form of a histogram. In Fig. 1, riji(s), which 
is the number of avalanches consisting of s grains for a given R, is presented for R = 3 in 
a semi-logarithmic plot. From the histograms of all the orifice sizes studied (1 < R < 4.5) 
it is possible to define the jamming probability J(R, N) as the probability of finding an 
avalanche of size less than or equal to N grains for a given R. That is to say, J(R, N) is the 
probability that an orifice of relative size R gets obstructed at least once before more than 
N grains fall through it. Then J can be evaluated as 

N 

E u r( s ) 

J(R,n) = s -^ (1) 

E n n( s ) 

We have obtained J for different values of R from the experiment; the results are shown 
in Fig. 2. The curves in the graph show that the bigger the size of the orifice, the smaller 
the probability that it gets blocked before an avalanche larger than N grains falls. Also, J 
grows with N for fixed R. These data are consistent with those found in a two dimensional 
hopper Q, obtained with a fixed nu mb e r of ^ains, m uch 8ma .,e r than in our expert. 

One could conceivably speculate from Fig. 2 that J would perhaps tend to a step function 
as iV — > oo. The question is whether for that limit the jamming probability is one below a 
certain size R c , while it is zero above R c . In other words, whether there is a finite value of R 
beyond which the orifice would not jam at all for an infinitely large silo. We can assess this 



3 



by evaluating the size of the orifice for which the jamming probability is one half -let us call 
it Rq(N). This can be easily obtained from the fit of J(R,N) (see Fig. 2). The particular 
fitting function for J chosen to obtain Rq(N) is almost irrelevant because the slope of J is 
quite large at R = Rq. In Fig. 3a we plot Rq versus N in semi- logarithmic scale. It is unclear 
whether Rq saturates. Several functions can be fit reasonably well to Rq(N) (see Fig. 3a), 
some of which saturate for N — ► oo, while others do not. Values of N much bigger than 10 5 
are experimentally inaccessible in our laboratory. Therefore we cannot conclude whether 
there is or not a critical value for Rq. The family of curves used to fit J also depends on 
a, which is related to the slope of J at Rq. The inverse of a is the width of the transition 
region from J = 1 to J = 0. The values for 1/a are plotted in Fig. 3b. It decreases with 
N but it is not possible to ascertain if it cancels out for a finite value of N, as would be the 
case for a step function. Both parameters, Ro and a, are not related in a simple way. 

These results can be understood in the framework of a one- dimensional percolation model. 
Let us call pr the probability that one grain gets past the outlet of size R without blocking 
it. Assuming that these events are independent, an avalanche of size s would consist of one 
jamming event, s non-jamming events and another jamming event. The probability of this 
sequence is 

n R (s) = (1 - PR ) 2 p s R (2) 

This also holds if one considers the grains falling in clusters, or groups, rather than one at 
a time. In this case, Pr is the probability that a group of grains passes the outlet without 
jamming, i.e. without forming an arch. If outgoing groups have k grains in average, the 
probability of finding an avalanche of s grains would be ur(s) = (1 — PR) 2 p S R k ' ■ From Eq. |2] 
one obtains logins) = 2 • log{l — Pr) + s ■ log(pR), so logijiR) has a linear dependence on s, 
exactly as observed experimentally. The value of Pr can be determined from the histogram. 
The exponential dependence of ur(s) on s precludes the interpretation of these data in 
the framework of self-organized criticality (SOC) because there exists a characteristic 
avalanche size manifested in the exponential distribution of events. Conversely, experiments 



with rough grains in different experimental conditions 



and numerical simulations in a 



two dimensional silo |3] yield a power law distribution - a signature of SOC. 

Incidentally, the exponential distribution is a compelling evidence that arch formation 
is an uncorrelated process. From the temporal point of view, this is also revealed by the 



autocorrelation function of the series of avalanches and by the first return map (see Fig. 1, 
inset), which show no sign of memory between consecutive avalanches. 

From the avalanche size distribution given by Eq. |2J the mean size avalanche S can be 
explicitly calculated, yielding: 
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o2 



Q _ s=0 _ 1 +PR / q \ 
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The mean avalanche sizes obtained from our data are shown along with the above ex- 
pression in Fig. 4. Indeed the histogram given by Eq . (J2J) and the mean avalanche size (JHJ) 
are those of a one- dimensional percolation model [14J, which leads to a second order phase 
transition, but only at pn = 1. This is marked by a divergence to infinity in the mean 
avalanche size S, which is a measure of the susceptibility of the system. Let us remark that 
the limit pr = 1 is unattainable, so such a phase transition would never be observed. In 
practice, however, pr could be made high enough (if the radius of the outlet orifice is big 
enough), so that jamming would seldom be observed. 

We can put forward the idea that the jamming transition bears an analogy with the glass 
transition. By increasing R one selects larger bridges only, so 1/R could be thought of as the 
"temperature". In the glass transition, the characteristic viscous time is certainly infinite 
for zero temperature; for a finite temperature, it can be extremely large but a power law 
divergence is not found In the same sense, the mean avalanche size (or the susceptibility, 
for that matter) is infinite for 1/R —> 0, and could be very large for finite values, but there 
would not exist a critical exponent for S. 

In this paper, we have shown that jamming events during the discharge of a silo of glass 
beads are uncorrelated and show no evidence of self-organized criticality, in contrast with 
other avalanche studies jl^l . Q] which were performed in different geometries. Moreover, 
we showed evidence that suggest that the classical claim of the existence of a critical outlet 
size beyond which no jamming occurs -i.e. J becomes a step function- could be invalid in 
a thermodynamic sense. Comparison between experiment and recent numerical simulations 
jl6 | could help to gather more understanding. Finally, other interesting aspects that have 
not been considered, such as the influence of the roughness of the beads and their packing 
fraction, also merit further investigation. 
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FIG. 1: Histogram for the number of grains that flow between two successive jammings. Data 
correspond to R = 3 (the beads have a diameter of 2 mm., and the circular orifice is 6 mm. wide). 
More than four thousand events have been recorded. The line is a linear fit. Inset: first return 
map, i. e. the avalanche size nt versus the next avalanche size rit+i- Note that t is just a correlative 
index ordering the sequence of avalanches. 
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FIG. 2: Jamming probability J as a function of R. Each series corresponds to a fixed TV, i.e. 
the probability of the orifice getting jammed before N grains pass through it. Data taken for 
several bead diameters are displayed. Succesive series, from left to right, correspond to N = 
3, 10, 30, 100, 300, 1000, 3000, 10000 and 30000 beads. The fit corresponds to a hyperbolic tangent 
with two adjustable parameters: J(R, N) = 1 ~ tanh [^( R R o)\ _ ^ wo parameters of this family of 
curves, a and Rq, are shown in Fig. 3. Error bars are statistical. 
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FIG. 3: (a) Rq versus N. Note the logarithmic scale on the horizontal axis, (b) 1/a as a function 
of N in semilogarithmic scale. Both Rq and a were obtained from the fits of Fig. 2. 
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FIG. 4: Mean avalanche size S as a function of pr\ experimental data (points) and model (solid 
line). For the experimental data, pr has been estimated as (jv^+aT) > where iVj, is the total number 
of fallen beads and Nj the total number of jammings, for a given orifice size. Inset: the same data, 
plot in logarithmic scale as a function of 1 — pr. 
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